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In these notes, we review some recent mathematical resuhs concerning the derivation of ef- 
fective evolution equations from many body quantum mechanics. In particular, we discuss the 
emergence of the Hartree equation in the so-called mean field regime (for example, for systems of 
gravitating bosons), and we show that the Gross-Pitaevskii equation approximates the dynam- 
1-^^ ' ics of initially trapped Bose-Einstein condensates. We explain how effective evolution equations 

, can be derived, on the one hand, by analyzing the so called BBGKY hierarchy, describing the 

time-evolution of reduced density matrices, and, on the other hand, by studying the dynamics of 
coherent initial states in a Fock-space representation of the many body system. 

^ ' 1 Introduction 

a^ 

0^ ■ Systems of interest in physics are typically characterized by a very large number of interacting degrees 
of freedom. Direct application of fundamental theories to study the properties of these systems is 
impossible. One of the main goals of statistical mechanics consists therefore in the derivation of 
simpler effective theories, providing a good approximation to the fundamental equations in the 
relevant regimes. 

We consider here a quantum mechanical system of particles moving in the three dimensional 
space and described by a wave function ip^ G L^(M^^,dxi . . . dx]\f). Physically, the absolute value 
squared \tpN{xi, ■ ■ ■ ,X]\i)\'^ determines the probability density for finding particle one at xi, particle 
two at X2 and so on. More generally, an observable of the system is a self-adjoint operator A acting 
on the L'^{M.^^,dxi ...dx^)', the physical expectation of this observable in the state described by 
the wave function V'Af is given by the L^-inner product {ipN, AtPn)- The positions of the particles 
are associated with multiplication operators (the operator A = Xj, acting by multiplication with 
Xj, measures the position of the j-th particle). The momenta of the particles are associated with 
diff'erential operators (the operator A = —iVx^ measures the momentum of the j-th particle). 

In nature, there exist two different types of particles, known as bosons and fermions, which behave 
differently with respect to permutations of the particles. The wave functions of bosonic systems are 
symmetric with respect to any permutation of the A'" particles (in the sense that tpNixni, ■ ■ ■ , x^n) = 
^N^xi, . . . ,xn) for any permutation vr G Sn). The wave functions of fermionic systems, on the 
other hand, are antisymmetric with respect to permutations (^Ar(x7ri, . . . , x-j^n) = '^■k'4'n{xi-, ■ ■ ■ ,xn), 
where a-j^ is the sign of the permutation vr G Sn). In these notes, we will exclusively consider bosonic 
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systems, and therefore we will always assume the wave function ■0jv to be symmetric with respect 
to permutations. We will denote by L^(M^^) the subspace of L^(M^^) consisting of all permutation 
symmetric functions. 

The evolution of the quantum system is governed by the A^-particle Schrodinger equation 

idt'ipN,t = HN'ipN,t (1-1) 

for the time dependent wave function ipN^t (the subscript t just indicates the time dependence of 
the wave function; time derivatives will always be written as dt). On the r.h.s. of the Schrodinger 
equation, Hj^ is a self-adjoint operator acting on Lg(M^^) and commonly known as the Hamilton 
operator (or Hamiltonian) of the system. It typically has the form 

N N 

H = Y, (-^-. + ^ext(x,)) + A^^ y(xi - Xj) (1.2) 

j=l i<j 

where T4xt is an external potential, A G R is a coupling constant and V{xi — xj) describes the 
interaction between particles i and j. Here, and in the following, we restrict our attention to Hamilton 
operators with two-body potentials, neglecting interactions which depend at the same time on three 
or more particles. 

The Schrodinger equation (jl.ip is a linear equation and its unique solution can be obtained by 
applying the unitary group generated by the self-adjoint operator H]\f to the initial wave function 
'4>N,t=0i i-e. V'Af,* = e~*^VAf,o- in other words, the well-posedness of the A^-particle Schrodinger 
equation p.ip is not an issue (although, depending on the properties of the potentials 14x1, V , it 
may not be easy to show that the unbounded operator H, defined on an appropriate domain, is 
self-adjoint). When studying systems of interest in physics, however, it is generally very difficult to 
extract qualitative or quantitative information (beyond its well-posedness) from (jl.ip . The reason 
is that, in typical situations, the number of particles involved in the dynamics is huge. It varies 
between N ~ 10^ in extremely dilute samples of Bose-Einstein condensates, to ~ 10^'^ in chemical 
samples, up to > 10^^ in stars and other astronomical and cosmological systems. For such values 
of A^, it is impossible to solve (|l.ip directly. Instead, it is desirable to derive effective equations 
approximating the dynamics governed by (jl.ip in the relevant interesting regimes. 

Probably the simplest non-trivial situation where it is possible to obtain an effective approxima- 
tion of the dynamics (jl.ip is the so-called mean field limit. The mean field regime is characterized 
by a large number of particles which interact very weakly with all other particles in the system, 
so that, effectively, the two-body interaction experienced by each particle can be approximated by 
an averaged, mean field, potential. The Hamiltonian (II. 2p describes a mean field regime if TV ^ 1 
(large number of particles) and A ^ 1 (weak coupling), so that A^A =: k remains of order one (this 
guarantees that the total potential, the total force acting on each particle is of order one). To study 
the evolution in the mean field regime we consider a factorized initial wave function = 
(meaning that '4)n{xi, ■ ■ ■ ,xn) = YljLi '■Pi^j)) for an arbitrary one-particle orbital G L^(]R^) and 
we study its evolution ij)N,t = e~^^'^^'ipN,o, as generated by the mean field Hamilton operator 

N 

Hn = Y^ (-^-. + VeA^j)) +^Y1 Vix, - X,) (1.3) 

j=l i<j 

in the limit of large A^. Because of the interaction, factorization cannot be preserved by the time- 
evolution. Nevertheless, because of the mean field character of the Hamiltonian, we may expect that 
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factorization is approximately preserved in the limit of large A^. If this is true, we would have 

N 

^nA^i^-'-^^n) ~Y\.y^t{xj) (1.4) 
i=i 

in an appropriate sense and for an appropriately evolved one-particle orbital ft G -L^(R^). Under 
this assumption, it is simple to find a self-consistent equation for the evolution of the orbital (ft- In 
fact, if ()1.4p holds true, the particles are approximately distributed in space independently of each 
others, with the probability density Iv'tP- The total potential experienced by, say, particle j, can 
therefore be approximated by 

We conclude that the one-particle orbital ipt must evolve according to the self-consistent, non- linear, 
Hartree equation 

idtipt = -A(pt + Vc^t^pt + i^(y * \'^t\'^)v^t ■ (1-5) 

To understand in which sense (11.40 can be correct, we introduce the concept of reduced density 
matrices associated with a bosonic wave function il^N,t- We define, first of all, the density matrix 
associated with Va^,* as the orthogonal projection 7jv^( = \^N,t){'>pN,t\ onto ipiy,t (we use here the 
shorthand notation |V')(V'I fo denote the orthogonal projection onto ip). For A; = 1, . . . , A^, we define 

(k) 

then the fc-particle reduced density matrix ^j^'^ by taking the partial trace of ^N^t over the last N — k 
variables, i.e. 

^nI = Trfc+i,...,iv 7N,t ■ (1-6) 

(k) 

In other words, the /c-particle reduced density matrix 7jv fc is defined as the non-negative, trace-class 



operator on L^(M^'^) with the integral kernel 

(A;) / _ / / \ 

77V t V'^^l ) • • • ) ^fc j 3^]^ , . . . Xf^) 

(1-7) 

dXk+l . ..dxN '4^N,t{xi, . . .,Xk,Xk+i,. . . ,XN)ll^N,t{x'i, ■ ■ . . . . ,Xiv) . 

(A;) (k) 

By definition Ti^j^^ = 1 for all k = 1,. . . ,N. Using the fc-particle reduced density 7]^^ one can 
compute the expectation of A;-particle observables (observables depending non-trivially on at most k 
particles). In fact, 

(^iv„ (OW l(^-^)) ^Tv,) = Tr (O(^) l(^-'=)) 7^, = Tr 0(^)7^ . 

It turns out that reduced density matrices provide the correct language to understand the con- 
vergence (|1.4|) in the limit of large N. 

Theorem 1.1. Under appropriate assumptions on the potentials V,Vext (see discussion below), 
consider the factorized N -particle wave function ipj\j- = tp®^ , for an arbitrary tp G H^{^). Let 
ipN,t = e~*^'^*VAf the evolution of ipN, ols generated by the mean field Hamiltonian lll.3\) . Then, 
for any fixed k £ N and t gW, we have, as N ^ 00, 

TSfl^i^iXv^^r' (1-8) 

where ipt is the solution of the nonlinear Hartree equation ( [i.5|) with initial data (pt=o = (p- The 
convergence in \1.8\) is in the trace class topology. 
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In particular, Theorem 11.11 implies that, for any bounded /c-particle observable 0^^\ and for any 
fixed t G M, we have 

as — 7- oo. In other words, as long as we compute the expectation of observables depending only on 
a fixed number of particles, the solution of the full Schrodinger equation ^l^N,t can be approximated, 
as N ^ oo, by the product ipf^ of solutions of the Hartree equation p.5|) . 

The first proof of this result was obtained by Spohn in [28], for bounded interaction potential 
V G L°°(M^). Later, Spohn's approach was extended by Erdos and Yau in [llj to interactions with 
Coulomb singularity V{x) = ±l/|x| (partial results for the Coulomb case were also obtained by 
Bardos-Golse-Mauser in [ij). Note that, in all these works, the important assumptions concern the 
interaction potential; only very minor conditions have to be imposed on the external potential Vext 
(to make sure that the operator ffjv is self-adjoint). 

The main idea of the approach developed in [28j and then extended in |T] is to study directly 
the time-evolution of the reduced density matrices. From the Schrodinger equation (II. ip for the 
wave- function il^N,t it is simple to derive a hierarchy, commonly known as the BBGKY hierarchy 
(BBGKY stands for Bogoliubov, Born, Green, Kirkwood, Yvon), consisting of N coupled equations 
describing the evolution of the reduced densities. For A: = 1, . . . , A^, we find 



+ 



Kj 



(1.^ 



+ ( ^ ~ ) ^ '^^^^^ ^^^^ ~ ^^+1^' W 



(fc+i) 



Here, we use the convention that 7^^"*"^^ = 0. In the last term, Trj^.^.! indicates the partial trace over 
the (A: -|- l)-th particle. It is simple to understand the origin of the three terms on the r.h.s. of (jl.Op . 
The first term corresponds to the kinetic energy of the first k particles. The second term describes the 
interactions among the first k particles. The last term, on the other hand, describes the interaction 
between the first k particles and the other (N — k) particles, which are integrated out in the definition 



of 7^^^ (this is why the last term depends on 7]y7^'' and not only on "y)^[)- As N ^ oo, the BBGKY 



(k+i) 



hierarchy converges, at least formally, to an infinite hierarchy of equations. Suppose that {7^\}fc>i 

denotes a limit point of the sequence {7^|}^i (with respect to an appropriate product topology). 
Then, from (jl.9p . we may expect the limit point to satisfy the infinite hierarchy (k > 1) 



A.^. + ^ext {xj ) , I + Tr^+i V{xj 



Xk- 



-l)>7oo,t 



(1.10) 



It is then worth noticing that this infinite hierarchy has a factorized solution. In fact, it is simple to 



check that the ansatz 7, 



(k) 

oo,t 



\ipt) {iptl^'' solves (jl.lOp if and only if (ft solves the nonlinear Hartree 



equation (jl.Sp . This observation suggests that, in order to obtain a rigorous proof of Theorem ll.H 
one can proceed as follows. First, one has to show the compactness of the sequence {7^{}^]^ with 
respect to an appropriate (weak) topology. In the second step, one has to identify the limit points 
{7oo\}a:>i of the sequence {7^{}^i as solutions of the infinite hierarchy (jl.lOp . Finally, one has to 
show the uniqueness of the solution of the infinite hierarchy (jl.lOp . These three steps immediately 
imply the convergence of the sequence {7^^}^! towards the factorized solution {\(pt) {^t\'^^} k>i with 
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respect to the weak topology; since \Lpt){ipt\ is a rank one projection, it is then easy to check that 
the weak convergence also implies strong convergence (in the trace- norm topology). Let us remark 
here that, as noticed by Chen and Pavlovic in [5], this approach can also be extended to many-body 
systems interacting via three-body potentials (namely potentials depending at the same time on 
three particles). 



2 Dynamics of Bose-Einstein Condensates 

A trapped Bose gas can be described by the Hamilton operator 



TV TV 

hT = E (-^-. + ^-t(x,)) + iVV(iV(x, - X,)) (2.1) 



acting on the Hilbert space L^(M^^, dxi, . . .dx^). Here V^xt is a confining potential (with Veyit{x) — ?■ 
oo as |x| — )■ oo), and y is a short range (possibly compactly supported), repulsive (meaning that 
V{x) > 0, for all X S M^), bounded interaction potential. In (j2.1|) . the interaction potential scales 
with the number of particles A^, so that its scattering length is of the order The scattering 

length is a physical quantity which measures the effective range of the potential V . It is defined 
through the solution of the zero-energy scattering equation 

-A + ^y(x))/ = (2.2) 

with the boundary condition f{x) — t- 1 as |x| — )• oo. For large f{x) has the form 

/(x)~l- -^ + 0(1x1-^) 
\x\ 

for an appropriate constant aoi which is defined to be the scattering length of V . Equivalently, the 
scattering length oq can be defined through 

Svrao = j dxV{x)f{x) . (2.3) 

It is simple to check that, for repulsive potentials, ao > 0. By scaling, it is also clear that, if 
Oo denotes the scattering length of the potential V , the scattering length of the rescaled potential 
yN{x) = N'^V{Nx) appearing in (j2.ip is given by a = a^/N . 

In [22], Lieb, Seiringer and Yngvason proved that, in the limit of large A, the ground state energy 
of the Hamiltonian (j2.ip can be approximated, for large A, minimizing the so called Gross-Pitaevskii 
energy functional. More precisely, they showed that, if Ep^ denotes the bottom of the spectrum of 
the operator (j2.ip . then 

lim = min SQ-pi'^) 

N^oo A i/jGL2{R3):||<^||=1 

where we defined the Gross-Pitaevskii energy functional 

egA^) = j dx (|V(^(x)|2 + yext(2;)|vp(x)p + ATia^\^{x)\^) (2.4) 

for arbitrary ip G i/^(IR^). In particular, in first order, the ground state energy depends on the 
interaction potential V only through its scattering length oq (the precise profile of the potential is 
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irrelevant). In [21], Lieb and Seiringer also showed that the ground state of (j2.ip exhibits complete 
Bose-Einstein condensation, in the sense that, if denotes the one-particle reduced density matrix 
associated with the ground state vector tpjsi of (|2.1|) . then 

tJ)^ ^ |0GP)(0GP| , asiV^oo, (2.5) 

where c^gp £ -L^(IR'^) denotes the minimizer of the Gross-Pitaevskii energy functional (12. 4p . Physi- 
cally, ()2.5p tells us that, in the ground state of ()2.ip . almost all particles (all particles up to a fraction 
vanishing as — )• 00) are in the one-particle state described by (pcP- 

The analysis of \22\ [2T] implies that if we prepare a trapped Bose gas at sufficiently small tem- 
peratures, the system will condensate into the minimizer of the Gross-Pitaevskii energy functional. 
It seems then natural to ask what happens if we perturb the gas, for example by removing the 
external traps (typically consisting of strong magnetic fields). The system will immediately react 
to the perturbation and it will begin to evolve. The evolution will be generated by the translation 
invariant Hamiltonian 

N N 

Hn = Y1 - a, . + N^ViN{xi - X,)) . (2.6) 

j=l i<j 



While the results of [221 [21] show the validity of the Gross-Pitaevskii theory for predicting the ground 
state properties of trapped Bose gases, the next theorem shows that it also correctly describes the 
evolution of initially trapped Bose-Einstein condensates. 

Theorem 2.1. Suppose V >0, V{-x) = V{x), \V{x)\ < (1 + \x\^y , for some fj > 5/2 and define 
H]\[ as in 112. 6\) . Consider a sequence of N -particle wave functions ipN G L^(M^^) such that 



• ijjN has finite energy per particles; there exists a constant C > with {ipM^Hi^ipisi) < CN . 

• ipN exhibits complete condensation; the one-particle reduced density matrix associated with 
ipN is so that 7^^ — )• for some ip G ff^(M^). 

Let il^N^t = e~*^'^*'i/'Af ■ Then, for every fixed A; G N and t G M, we have 



Tr 







as N ^ 00, where (ft is the solution of the time- dependent nonlinear Gross-Pitaevskii equation 

idt^t = -^Vt + 8irao\<pt\'^ipt (2.7) 

with initial data (pt=o = '■P- Here oq denotes the scattering length of the (unsealed) interaction 
potential V. 

This theorem, which was proven in El O [10] , states that complete condensation is preserved by 
the time- evolution and that the condensate wave function evolves according to the Gross-Pitaevskii 
equation ()2.7p . Note that a different derivation of the Gross-Pitaevskii equation was proposed by 
Pickl in [261. 

The Hamilton operator (j2.6p can be written as 



N N 



N 

j = l Kj 
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and one may try to interpret Hi\i as a mean field Hamiltonian, with interaction potential fAr(x) = 
N^V{Nx). This observation suggests that, in order to show Theorem 12.11 one can use again the 
strategy outlined at the end of Section [H It turns out, however, that one should be very careful with 
this analogy, and that, although the general strategy based on the study of the BBGKY hierarchy 
still applies, several changes are needed and here a much deeper analysis of the A^-particle dynamics 
is required. 

Let us briefly explain the new challenges appearing in the proof of Theorem 12.11 as compared 
with the mean field regime discussed in Section [TJ Since, formally, vp^{x) — )• bo6{x), with Bq = 
J dxV{x), the naive analogy with the mean field situation suggests that the dynamics generated by 
the Hamiltonian (j2.6p can be approximated, for large N, by the Hartree equation 



This equation has the same form as the Gross-Pitaevskii equation ()2.7p . but the wrong coupling 
constant in front of the non-linearity. From the physical point of view, it is not surprising that the 
mean field analogy fails. The mean field regime is characterized by a large number of very weak 
collisions among the particles. The dynamics generated by (12. 6p . on the other hand, is characterized 
by rare and, at the same time, very strong collisions (particles only interact when they are at 
distances of order 1/A^ from each others). As a consequence, it turns out that the correlations 
among the particles, which were negligible in the mean field regime, play here an important role 
and, in particular, are crucial to understand the emergence of the scattering length oq in the Gross- 
Pitaevskii equation (j2.7p . 

Because of the singularity of the interaction potential, the solution of the Schrodinger equation 
'ipN,t = e~'^^'^^-4>N develops a short-scale correlation structure which lives on the same length scale 
1/A^ characterizing the potential. The singular correlation structure can be approximately described 
by the zero-energy scattering equation f{x) defined in (12. 2p (more precisely, the singular structure 
is described by the solution fN{x) = f{Nx) of the zero-energy scattering equation with rescaled 
potential N'^V{Nx)). So, the effective, average potential which is experienced say, by particle j due 
to the interaction with the other (A^ — 1) particles can be approximated here by 



where we used the characterization (j2.3p of the scattering length oq. Hence, taking into account 
the correlations, we obtain the Gross-Pitaevskii equation ()2.7p with the correct coupling constant 
in front of the nonlinearity. This heuristic argument shows that one of the main new challenges 
consists in proving that the A^-body dynamics really develops a singular correlation structure which 
can be described, in good approximation, by the solution of the zero-energy scattering equation (j2.2p . 
Another major challenge, compared with the results obtained in the mean field regime, is the proof of 
the uniqueness of the infinite hierarchy. The main problem here is that the interaction potential given 
in the limiting hierarchy by a delta-function cannot be controlled, in three dimensions, by the kinetic 
energy. As a consequence, uniqueness is proven in [7J by expanding the solution in a complicated 
diagrammatic expansion in terms of Feynman graphs; to control the many contributions in this 
expansion, it is very important to use the smoothing effects of the free evolution, which effectively 
regularize the singular interaction potential. A new and simpler proof of the uniqueness of the infinite 
hierarchy was obtained by Klainermann and Machedon in [18] (later, this approach was extended 
by Chen and Pavlovic in These works show the uniqueness of the infinite hierarchy in a class 



idtift = -A(pt + {bo6 * {iptl"^) ipt = -A(pt + bo\'^t\'^'^t ■ 
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of densities satisfying certain space-time bounds. Unfortunately, so far it has not been possible to 
show that the limit points of the sequences {'y^j^\}^^i satisfy these bounds; as a consequence, so far 
it was not possible to apply the results of |18l dj to prove Theorem 12.11 (in one and two dimensions, 
on the other hand, the results of [TBI S] can be applied to show the analogous of Theorem 12. Ij see 
p!7] ) . For more details of the proof of Theorem 12. 1^ we refer to [3 [U [9l [TO] . 



3 Mean Field Evolution of Coherent States 

The main drawback of the methods outlined in Sections [1] and [2] is the lack of precise bounds on the 
difference between the many body dynamics and the effective Hartree evolution. With an expansion 
of the solution of the BBGKY hierarchy (II. 9p . it is possible to show (at least for bounded interaction 
potentials) that there exists constants C, Tq > such that 



Tr 



l^^l-\v>t){^<^ (3.1) 



for all t G R with \t\ < Tq. Similar bounds can also be obtained for the reduced fe-particle densities, 
for fixed k > 2. Unfortunately, (|3.1|) is only valid for short times; for t > Tq, one can still iterate 
(j3.1|) , but one only obtains much weaker estimates of the form 



Tr 



It turns out that one can derive stronger (optimal) bounds on the rate of the convergence of 
the many body evolution towards the mean field Hartree dynamics using techniques originating in 
quantum field theory. These techniques, which were first introduced by Hepp in [16] for the analysis 
of the classical limit of quantum mechanics, and later extended by Ginibre-Velo in J13J, are based 
on a Fock-space representation of the bosonic systems and on the study of the dynamics of so-called 
coherent states. A different approach, which also leads to explicit bounds on the rate of convergence 
was proposed by Pickl in [25] , and then applied by Knowles and Pickl to the derivation of the Hartree 
equation with singular potentials in [19]. 

The bosonic Fock-space over L^(M'^) is defined as the direct sum 



jr = c e Ll{W^'', dxi,... dxn) 



n>l 



where Lg(M^"') denotes the subspace of L^(R^") consisting of functions symmetric with respect to 
permutation of the n particles. Vectors in the Fock-space are sequences „>o, where 

is an n-particle bosonic wave function. The idea behind the introduction of the 
Fock-space is that we want to study states where the number of particles is not fixed. Clearly, T has 
the structure of a Hilbert space with the inner product 



n>l 



The vector Q. = {1, 0, 0, ... } G is called the vacuum and describes a system with no particles. An 
important operator on T is the number of particle operator AA, which is defined by 
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The vacuum Q is an eigenvector of with eigenvalue zero. More generaUy, vectors of the form 
{0, . . . , 0, -0^™), 0, . . . }, having a fixed number of particles, are eigenvectors of (with eigenvalue m). 
On T, we define the Hamilton operator by 'HN{'ip^"'^}n>i = {'HN^'^^"'^n>i, with 



j = l i<j 

By definition, T-L]\f leaves each n-particle sector J^n (defined as the eigenspace of M associated with 
the eigenvalue n) invariant. Moreover, on the A^-particle sector, T-L]\f agrees with the mean field 
Hamiltonian Hj\f defined in In particular, if we consider the Fock-space evolution of an initial 

vector with exactly particles, we find 

e-^«^*{0, . . . , 0, Viv, 0, . . . } = {0, . . . , 0, e-^^^ViV, 0, . . . } 

exactly as in Section [TJ The advantage of working in the Fock space is that we have more freedom in 
the choice of the initial data. We will use this freedom by considering a class of initial data, known 
as coherent states, with non-fixed number of particles. 

It is very useful to introduce, on the Fock space creation and annihilation operators. For 
/ G L'^{M.^), we define the creation operator a*{f) and the annihilation operator a(/) by setting 

1 " 

(«*(/)V')^"^ {Xl, . . . ,Xn) = —^Y] fixj)'^^'^'''^Hxi, Xj^i,Xj+i, ...,Xn) 
(a(/)V^)(") (Xi, . . . , X„) = [ dx Ji^i^^^+^Hx, XI,..., Xn). 



The operators a*{f) and a{f) are densily defined and closed. It is easy to check that, as the notation 
suggests, a*{f) is the adjoint of a{f). Creation and annihilation operators satisfy the canonical 
commutation relations 

[a{f),a*{g)] = {f,g) and [a{f),a{g)] = [a* {f) , a* (g)] = 

for any f,g & L^(]R^) (here {f,g) denotes the L^-inner product). Physically, the operator a*{f) 
creates a particle with wave function /, while a(/) annihilates it. As a consequence, a state with N 
particles all in the one-particle state if, can be written as 

{o....,o.^«o....} = <«n. 

It is also useful to introduce operator- valued distributions ax, a* defined so that 



a{f) = J dxf{x)ax, and a*{f) = J dxf{x)a*. 

With this notation, a* ax gives the density of particles close to x € M^. The number of particles 
operator can formally be written as 

N = 1 dx a„.ax ■ 



Similarly, the Hamilton operator H.^ can be formally expressed in terms of operator- valued distri- 
butions as 

Hn = J C^xa*x Vxax + V;xt(x)a*ax-) + ^ J ^^'^V ^(^ ~ y)(i*xa*yayax ■ 
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The fact that every term in the Hamiltonian contains the same number of creation and annihilation 
operators means that T-Lm commutes with the number of particles or, equivalently, that the number 
of particles is preserved by the time-evolution enerated hy T-Ln. 

For later use, we observe that the creation and annihilation operators are not bounded; however, 
they can be bounded by the square root of the number of particles operator, in the sense that 

\Hm\ < 11/11 IIAAi/'V'll and ||a*(/)V^|| < ||/|| ||(AA + l)i/2^|| (3.2) 

for every ip £ J-, f G L^(M^) (here ||/|| indicates the L^-norm of /). 

As mentioned above, we are going to study the evolution of initial coherent states. For arbitrary 
If G L^(R^), we define the Weyl operator 

Wif) = e»*(^)-"(v) . 

The coherent state with wave function ip is then defined as W{(p)i}. The Weyl operator W{(p) is a 
unitary operator; therefore the coherent state W{(p)Q always has norm one. Moreover, since 



j=Q 

the coherent state W{(p)Q does not have a fixed number of particles. One can nevertheless compute 
the expectation of the number of particles in the state W{ip)il,; it is given by 

{w{ip)Q,Mw{ip)n) = yf. 

More precisely, it turns out that the number of particle in the coherent state W{ip)il. is a Poisson 
random variable with expectation and variance Hv^lp. The main reason why coherent states have 
nice algebraic properties (which will be used later on in the analysis of their evolution) is the fact 
that they are eigenvectors of all annihilation operators. Indeed 

a{f)W{ip)n = {f,ip)W{ip)^} 

for every f,(pG L^(R^). This is a simple consequence of the fact that Weyl operators generate shifts 
of creation and annihilation operators, in the sense that 

W*{ip)a{f)Wiip) = aif) + if, if), and W*{ip)a*if)Wiip) = a*{f) + {ip, /) . (3.3) 

Next, we study the evolution of an initial coherent state with respect to the dynamics generated 
by the Hamilton operator T-Ln- To reproduce the mean field regime analyzed in Section [H we choose 
the initial coherent state to have expected number of particles equal to N (the number of particles 
cannot be fixed, but at least we should fix its average to be given by N; otherwise the resulting 
evolution will not have anything to do with the mean field Hartree dynamics). 

Theorem 3.1. Suppose that the interaction potential V is such that, as an operator inequality 
< {1 — A). For (f G H^{W^), consider the initial coherent state 



Let "^1^,1 = e~'^'^^W{VNip)Q, and let T)^\ denote the one-particle reduced density associated with 
^Nt- Then there exist constants C,D >0 such that 



Tr 



r^^l-\ipt){ft\ <C^ (3.4) 
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for all t G M. Here (ft denotes the solution of the nonlinear one-particle Hartree equation l[1.5\) with 
initial data ft=o = f- 

Observe that the operator inequahty V'^ < C(l — A), which means that 

dxV{x)\ip{x)\'^ < C\\ip\\Hi 

for ah if G H^{W^), is satisfied for potentials with Coulomb singularities V{x) ~ ibl/|x|. 

From the convergence towards the Hartree dynamics for the evolution of initial coherent states, 
one can deduce a similar result for the evolution of initially factorized states with a fixed number of 
particles. To this end, one can use the fact that 

{0, . . . , 0, <^^^, 0, . . . } = ^ n = dNPNW{VNip)n 

where Pn is the orthogonal projection onto the A^-particle sector of and dN = e^^'^N~'^^'^\fN^. ~ 
N^/'^. Alternatively, one can write 

{0, . . . , 0, v?®^, 0, . . . } = / de e"*^ V(\/iVe^V)f^ 
to express the factorized state as a linear combination of coherent states. 

Corollary 3.2. Let the potential V be so that V'^ < {1 — A) and define the Hamiltonian as in 
U.3\) . Let -ipN = V'®^ o.'i^d ipN,t = e~'^^'^^ipN ■ Then there exist C,D > such that 

Tr lZ-\^t){^t\ 

for all t gM. Here <pt is the solution of the Hartree equation U.5\) . with initial data (pt=o = ^■ 

The details of how Corollary 13.21 follow from Theorem 13.11 can be found in |27l [2]. In the following 
let us briefly present the main ideas behind the proof of Theorem 13. li The main idea is to use the 
fact that the evolution of a coherent state remains approximately coherent. As we will see, it is 
possible to extract the coherent part of the evolved state, and then to focus on the evolution of 
the fluctuations, which, thank to the algebraic properties of the coherent states, can be expressed 
in simple and compact form. Let us remark that, in [141 I15j . Grillakis, Machedon and Margetis 
proposed to approximate the evolution of the coherent state not just by a coherent state but by a 
larger manifold of so-called Bogoliubov states; this allows them to obtain more precise approximation 
of the many-body evolution (this approach has been extended to systems with three-body interactions 
by Chen in [3]). 

The first observation is that the one-particle density matrix r^-*^ associated with the Fock-space 
state ^N,t has the integral kernel 

^N^i^'^y) = t;f, Tr^, — T {'i>N,t,alay^N,t) = ^{'i' N,t, a^ay^ N,t) ■ 

Expanding a* and Oy around their mean field values ^fNipt{x)^ VNiptiu), we obtain 

= ^ (n, W*{VN^)e'^^' - \/iV(^t(x)) (ay - VN^tiv)) e-'^^'W{VN^)n 
+ (Q,W*iVNip)e'^^^ (a* - VNip^{x)) e-'^^^W{^/Nip)Q 



11 



Using (j3.3p . we rewrite the last equation as 



\/ N \ 

+ (n,W*{VNip)e''^^^W{VNipt)alW%VNipt)e-'^^^ . 



Introducing the two-parameter group of unitary transformations 

U(t; s) = VF*(\/iV(/7t)e~^^'^(*-")VF(\/iV(/7,), (3.5) 

we find 



(3.6) 



N 



Let us first consider the first term on the r.h.s. of the last equation; recalling the bounds (j3.2p . we 
conclude that the contribution of this term to the I.h.s. of (|3.4p can be controlled by the r.h.s. of 
p.4p if we can control the growth of the expectation of the number of particles operator J\f with 
respect to the fluctuation dynamics U{t; 0), i.e. if we can prove that 

{n,u*{t;o)Aru{t;o)n) < Ce^i*i . (3.7) 

It is worth noticing that the fluctuation dynamics V({t; s) satisfies the Schrodinger equation 

idtU{t; s) = C{t)U{t; s), with U{s; s) = I 

with the generator 

C{t) = j dx [V xalV + Vey±{x) ala^, + {V * \'ft\^){x)ala:^) + j dxdyV{x-y)ipt{x)Tpt{y)alay 
+ ^ I dxdyV{x-y) {(ft{x)(pt{y)ala*y+lft{x)lpt{y)axay) 



1 



+ — / dxdyV{x - y) a* {aliptiv) + ayiptiv)) 



+ ^ y dxdy V{x - y) alalaya^ . 

(3.8) 

In contrast with the Hamiltonian Hn, the generator C{t) contains terms (the terms on the second 
and third line) which do not commute with the number of particles operator J\f (because the number 
of creation operators does not match the number of annihilation operators). Hence, not surprisingly, 
the expectation of M is not preserved by the fluctuation dynamics U{t;0). Nevertheless, it turns 
out that, if the condition V'^ < C{1 — A) is satisfied, it is still possible to control the growth of the 
expectation of Af, and to prove the bound (|3.7p . For the details, we refer to [27j . 
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The analysis of the second and the third term on the r.h.s. of (|3.6p requires shghtly more works. 
Formahy these terms seem to be of the order A^~^/^. To show that they really are of the order N~^, 
and prove the bound (j3.4|) . one needs to compare the fluctuation dynamics U{t;0) with a limiting 
dynamics Uoo{t;0), satisfying the Schrodinger equation 

idtUooit; s) = £oo(i)^oo(t; s) 

with the generator 

^oo(t) = j dx iy ^alS/ xax + Ve^t{x) ala^ + {V * \ipt\^){x)alax) + j dxdyV{x - y) ipt{x)Tpt{y)alay 

+ dxdyV{x-y) {ipt{x)ipt{y)ala*y +Tpt{x)Tp^{y)axay) 

obtained by formally letting — )• oo in (j3.8p . The point is that if we replace ly({t; 0) by Uoo{t] 0) in 
the second and in the third term on the r.h.s. of (j3.6p . these terms vanish. In fact 

{nMoo{t;Q)*axUoo{t;^)^) = {n,Uoo{t;OyalUoo{t;0)n)=0 

because, although Uodt; 0) does not preserves the number of particles, it does preserve the parity (in 
the sense that it commutes with the operator (—1)-^). This observation implies that, in the second 
and third terms on the r.h.s. of (|3.6p . the unitary evolution U{t; 0) can be replaced by the difference 
U{t; 0) —Uoo{t; 0); this produces the additional factor N^^/'^ (because the difference between the two 
generators, C{t) and Coo{t), is of this order), and explains why also the second and the third term 
on the r.h.s. of (j3.6p are of the order A^~^. Again, we refer to [2Tj for further details. 



4 Gravitational Collapse of Boson Stars 

As an application of the techniques discussed in Section [3l the last part of these notes is devoted 
to the study of the dynamics of so-called boson stars. These are systems of bosons with relativistic 
dispersion law interacting through classical Newtonian gravity (such particles are usually called 
semi-relativistic, or pseudo-relativistic bosons). The Hamilton operator has the form 

N N 

We are interested here in the mean field regime, characterized by ^ 1 and G ^ 1, so that k := NG 
remains fixed, of order one. Since, in the units we use, the gravitational constant G is approximately 
given by G ~ 10~^^, this means that we are interested in systems with approximately ~ 10^^ 
particles. To analyze this regime, we consider the time-evolution generated by 

N N 

if.=5:yrr^-^E^ (4.1, 

j = l Kj ■> 

in the limit N — )• oo. Inspired by the results discussed in the previous sections, we expect that, 
in this limit, the evolution generated by Hfq on factorized initial data can be approximated by the 
mean field semi-relativistic Hartree equation 

1 
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It turns out that the system under consideration is critical. This follows from the observation 
that, for large momenta, the kinetic energy \/l — A ~ |V| scales like the potential energy, as an 
inverse length. This implies that, for arbitrary N £ N, there exists a critical coupling Kcr{N) > 
such that 

inf {iPn,HniPn) > 

^/']VeL2(K3iV)J|^^|| = l 

for K < Kcr{N), and 

inf {ipN, Hn'^Pn) = -oo 



for K > Kcr- In other words, for small k > 0, the kinetic energy controls the potential energy, and 
the total energy is bounded below. For n > Kcr, on the other hand, the potential energy dominates 
and leads the total energy to arbitrary negative values. Criticality can also be observed on the level 
of the semi-relativistic Hartree equation (|4.2p . As proven in |20] . (j4.2p is locally well-posed in the 
energy space i7^/^(M^) for arbitrary coupling constants. Its global behavior, however, depends on 
the value of k. There exists namely a critical coupling Kcr > with the following properties. For 
K < Kcr, (|4.2|) is globally well-posed in H^^'^{Mj^), in the sense that all local solutions extend to all 
times. For k > on the other hand, ()4.2p has local solutions ipt G C{[0,T), H^^'^{M^)) exhibiting 
blow up in finite time, in the sense that 

||(/7t||j:^l/2 OO 

as t — 7- T~ . In [12], Frohlich and Lenzmann proved that arbitrary spherically symmetric initial data 
with negative energy lead, if k > Kcr to blow up in finite time (the spherical symmetry is believed to 
be just a technical condition). In the physics literature, the blow up solutions of the semi-relativistic 
Hartree equation (j4.2p have been used to describe the phenomenon of stellar or gravitational collapse, 
first predicted by Chandrasekhar. 

Mathematically, it seems important to understand whether the relation between the many body 
evolution and the semi-relativistic Hartree dynamics (14. 2p can be established rigorously. The first 
step in this direction was accomplished by Lieb and Yau; in [23], they proved that Kcr{N) — )• Hcr as 
N — )• OO, and that, for k < Hcr, the ground state energy per particle converges, as — )■ oo, to the 
minimum of the Hartree energy functional 



eHartree(v) ^ j ~ ~ ^ / ^^"^^ 



\x - y\ 

over all if £ L2(]R3) with \\ip\\ = 1. This proves that the Hartree theory successfully predicts the 
ground state properties of the boson star. Does it also describe their time-evolution? 

A first positive answer to this question was obtained in [6j, for the subcritical case k < k^. To 
consider the supercritical case k > k^, we first have to give a precise mathematical definition of the 
many body evolution. In fact, for n > Hcri^), it is not so simple to define the many body evolution 
because Hjsr is not bounded below. As a consequence, Hjsf does not have a unique and natural 
extension as a self-adjoint operator on the Hilbert space L^(M^^). Therefore, it is not clear how to 
define the unitary group e~^^^^ and it is not clear how to solve the many body Schrodinger equation. 
A possible way to avoid this problem is to consider weak solutions of the many body Schrodinger 
equation. Here, we follow a different approach. To circumvent the fact that the Hamiltonian is 
unbounded from below, we introduce an arbitrarily small, A'^-dependent cutoff a at in the interaction, 
defining the regularized Hamiltonian 

N 

HZ = y .Fi^^. - -Y , ^ — . 
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The cutoff Oat is assumed to be strictly positive for all G N, and to vanish in the limit of large 
A^. For arbitrary N £N, Hfj is now bounded below and generates therefore a well-defined unitary 
evolution on L?{^^). On the other hand, since a^v vanishes as N ^ oo, we do not expect it to 
considerably affect the macroscopic properties of the dynamics. From the physical point of view, the 
introduction of the A^-dependent cutoff is justifiable by the fact that, anyway, at very small distances, 
Newton's gravity is effectively regularized by the presence of other forces. 

Now, we are ready to study the relation between the many body evolution generated by H'^ and 
the semi-relativistic Hartree dynamics (j4.2p . The next theorem, proven in [23], shows that, if the 
nonlinear dynamics is well defined in a time interval [— T, T], then, in this time interval, (14. 2p really 
approximates the many body evolution in the limit of large N . 

Theorem 4.1. Let ip G H'^{R^), qat > such that ^ as N ^ oo, TpN,t = e'^^^^V^^- Let ipt 
be the solution of the semi-relativistic Hartree equation 1^4-^ ! with initial data ft=o = f- Fix T > 0, 
and assume that 

A := sup ||93t||Hi/2(R3) < oo . 

|i|<T 

Then, for every A; G N, there exists a constant C = C{k,T, X) > such that 

l^Nl-\Vt){^r' <c(j^^+a^) . (4.3) 



sup Tr 

\t\<T 



What happens now if the Hartree dynamics exhibits blow-up. The next theorem, which was also 
proven in [24J, shows that, if the Hartree dynamics blows up at time T, then also the many body 
evolution collapses if t — )• T, and, simultaneously, — )• oo (at a sufficient fast rate). By collapse of 
the solution 'ipN,t of the many body Schrodinger equation, we mean the following. The kinetic energy 
per particle, given by 

{i^N,t, (1 - A, Ji/Viv,t> = Tr (1 - A)i/27^^) , 
is finite, uniformly in A^, for alH < T, while it diverges to infinity if t — t- T and A^ — t- oo. 

Theorem 4.2. Assume that, in the definition of H'^, > N"^ for some (arbitrarily large) £ > 0. 
Let ip G H'^{M?), il^N,t = e~^^N^ip®^ . Let ipt be the solution of the semi-relativistic Hartree equation 
1^4-^ with initial data ipt=o = V- Assume that ipt is locally well posed in //^/^(M^) in the interval 
[0,T) and that it blows up at time T, that is \\<pt\\{ji/2 < oo, for all t G [0, T) and 

lim llc/^t II rri/2 = oo . 

Then, for every < t < T , we have 

2V(1-A)V2^« <oo (4.4) 
uniformly in N . Moreover, if N{t) — )• oo as t ^ T~ sufficiently fast, we have 

^hm Tr(l-A)V2^«^^^ = oo. (4.5) 

This theorem establishes the (dynamical) gravitational collapse of the boson star directly on the 
level of the many body evolution, justifying the use of the semi-relativistic Hartree equation (|4.2p 
for the description of the dynamics. 

The proof of Theorems 14. II and 14.21 relies on the ideas discussed in Section [3l In particular, we use 
a Fock-space representation of the system and we study the evolution of initial coherent states. The 
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main new challenge is that, in order to prove Theorem 14.21 (in particular, in order to show (j4.4p ). we 
establish the convergence of the one-particle reduced density associated with the evolution of the 
initial coherent state towards \(pt){^t\ in an energy norm (in Section O convergence was established 
in the trace norm, which, in the language of density matrices, is the equivalent of an L-^-norm). More 
precisely, we show that, for any < T (recall that T is the blow-up time of the nonlinear Hartree 
equation (02])), 

where the constant C > depends only on to and on sup^g[o,to] ll'/'tll_ffi/2- In Section [3l the proof 

of the convergence of the reduced density matrix F^''^ in the trace norm reduced to the problem of 
controlling the growth the number of particles operator with respect to the fluctuation dynamics 
U{t;0) defined in (|3.5p . Similarly, the proof of (|4.6p reduces to the problem of controlling the 
growth of the expectation of the (relativistic) kinetic energy operator with respect to the fluctuation 
dynamics. The fundamental reason why this is possible is that, after factoring out the (super-critical) 
Hartree dynamics from the (super-critical) many body evolution, the dynamics of the fluctuation is 
sub-critical. For further details, we refer to |24j . 



sup Tr 
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